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O ' Abstract 

o ■ 

' This paper concerns the reconstruction of the absorption and scattering 

Q_i! parameters in a time-dependent Hnear transport equation from full knowledge 

I of the albedo operator at the boundary of a bounded domain of interest. We 

present optimal stability results on the reconstruction of the absorption and 
scattering parameters for a given error in the measured albedo operator. 
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^ 1 Introduction 

■ Inverse transport theory has many applications in e.g. medical and geophysical 
tt^. imaging. It consists of reconstructing optical parameters in a domain of interest from 

! measurements of the transport solution at the boundary of that domain. The optical 
^ ' parameters are the total absorption (extinction) parameter a{x) and the scattering 
\Q parameter k{x,v',v), which measures the probability of a particle at position x to 
O ■ scatter from direction v' to direction v. 

Q The domain of interest is probed as follows. A known flux of particles enters 
the domain and the flux of outgoing particles is measured at the domain's boundary. 
O ' Several inverse theories may then be envisioned based on available data. The least 
^ favorable situation is when the density of outgoing particles is angularly averaged, 
' which means that only the spatial density of particles may be estimated and not the 
phase space (space and direction) density. Angular averaging may be necessitated 
by equipment cost, time of acquisition of the measurements, or low particle counts. 
5^ ] For uniqueness and stability results in this setting, we refer the reader e.g. to Bal 
and Jollivet [BJ2], Bal et al. [BLM], and Langmore [L]. 

A much more favorable situation is when the density of outgoing particles is 
angularly resolved. We may then be able to sample the outgoing distribution of 
particles as a function of time if sufficiently accurate equipment is available. In 
many setting however, only time independent measurements are feasible. 

The uniqueness of the reconstruction of the optical parameters from knowl- 
edge of angularly resolved measurements both in the time-dependent and time- 
independent settings was proved in ChouUi and Stefanov [CSl, CS2]. We also refer 
the reader to Stefanov [S] for a review of uniqueness results in inverse transport the- 
ory. Stability in the time-independent case has been analyzed in dimension d = 2,3 
under smallness assumptions for both optical parameters by Romanov [Rl, R2] and 
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in dimension d = 2 under smallness assumption for the scattering parameter by Ste- 
fanov and Uhlmann [SU]. Partial results on the stability of the reconstruction in the 
time-independent setting in dimension d = 3 were obtained in Wang [W] without 
smallness assumptions. Complete stability results in the time-independent case in 
dimension d > 3 were obtained by the authors in [BJl]. The present paper proves 
stability results for the time-dependent inverse transport problem. We restrict our- 
selves to the case of elastic scattering, where the velocity space may be modeled by 
the unit sphere S*^"^. Optimal results on the stability of the optical parameters are 
obtained in all dimensions d >2. 

The rest of the paper is structured as follows. Section [2] recalls useful results 
on the time- dependent linear transport equation. The main stability results of this 
paper are stated in section [HI They are based on a decomposition of the albedo 
operator used in [CSl] and recalled in section 13. 2[ Useful regularity results on the 
decomposition are stated in Proposition 3.2 and proved in section 4. Our first 
stability result is stated in Theorem 3.1. It shows how the Radon transform of 
the absorption parameter and a weighted norm of the scattering coefficient may 
be stably reconstructed from knowledge of the albedo operator. Under additional 
regularity assumptions, Theorem 3.2 shows the stability of the reconstruction of 
both optical parameters. Both stability results are proved in section 5. 

2 The forward problem 

In this section we introduce some notation and recall known facts about the well- 
posedness of the forward transport problem. 

2.1 The linear Boltzmann transport equation 

Let X be a bounded open subset of M'^, d > 2, with a boundary dX. We 
denote the diameter of X by diam(X) (diam(X) := sup(^^^y^^x^ 1^ ~ vD- Let i/{x) 
denote the outward normal unit vector to dX at a; G dX. Let F-i- = {{x,v) G 
dX X E>'^~^ I ± i'(x)v > 0}. For {x, v) E X x S'^^^ we define t±{x, v) and t(x, v) by 
T±{x, v) := inf{s G (0, +oo) | x ± sv ^ X} and r(x, v) := t_{x, v) + t^{x, v). 

Consider a : X x S'^'^ ^ M and A; : X x S'^-i x S'^-i ^ M two nonnegative 
measurable functions. We assume that (cr, k) is admissible when 



Let T > 77 > 0. We consider the following linear Boltzmann transport equation 
with boundary conditions 



— {t,x,v)+vVxu{t,x,v)+a{x,v)u{t,x,v)= / k{x,v',v)u{t,x,v')dv', {t,x,v) E {0,T)xXx'B'^ 



< a G L^{X X §^-1), 

< k{x,v',.) E L\S'^-^) a.e. {x,v') E X x 



(2.1) 





(2.2) 



u\{o,T)xr^{t,x,v) = (l){t,x,v), 
u{0,x,v) = 0, {x,v) EX xS' 
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where (t) G L^((0, T), L^(r_, d^)) and supp0 C [0,77]. 

We assume here that scattering is elastic, which imphes that the speed of the 
particles is preserved by scattering while only the direction of propagation may 
change. Elastic scattering is a good approximation in many applications in medical 
and geophysical imaging. Our results are stated for a (normalized) velocity space 
equal to the unit sphere S'^"^. Generalizations to other velocity spaces may be 
obtained as in e.g. [BJ2] and [CSl, CS2]. 

2.2 Semigroups and unbounded operators 

We introduce the following space 

Z := {/ e L\X X §<^-i) I vVJ e L\X X E."-')}, (2.3) 
WfWz := ||/|Ui(XxSd-i) + lbV^/||Li{XxSrf-i); (2.4) 

where vVx is understood in the distributional sense. 

It is known [CI, C2] that the trace map /3_ from C\X x S'^-^) to C{T_) defined 

by 

/?-(/) = /|r_ (2.5) 

extends to a continuous operator from Z onto L-'^(r_, t+(x, v)d$,{x, v)) and admits a 
continuous lifting. Note that L^{T_, d^) is asubset of the space L^(r_, t+(x, v)d^{x^ v)). 

We introduce the following notation 

AJ = -af, A2f = / k{x, v', v)f{x, v')dv'. (2.6) 

As ((T, k) is admissible, the operators Ai and A2 are bounded operators in L^{X x 
S'^-i). 

Consider the following unbounded operators 

TJ = -vVJ + A^f, D{T^) = {/ e Z I /|r_ = 0}, (2.7) 
Tf = Ti/ + A^f, D{T) = L>(Ti). (2.8) 

The unbounded operators Ti and T are generators of strongly continuous semigroups 
Ui{t) and U{t), respectively, in V-{X x S*^"^) (see e.g. [DL, Proposition 2 p.226]). 
In addition, Ui{t) and U{t) preserve the cone of positive functions and Ui{t) is given 
explicitly by the following formula 

Ui{t)f = e-fo^^^-'^'^^dsjf^^ _ i^^y^Q^r^ _ fy^^^^ foj- a.e. {x,v) eXx (2.9) 
for f eL\X X S'^-i), where 

e{x,y) = (^'^'' + - ^ .^"^ P ^ (2.10) 
^ ' \ otherwise, ^ ^ 

for (x, y) eW^ X W^. 

We will use the Duhamel formula 

U{r')^Ui{r')+ Ui{r' - s')A2U{s')ds', ior r' >0. (2.11) 
Jo 
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2.3 Trace results 

We introduce the following space 



W := |m G L\{0,T) X X X S^^-^) | + u G L\{0,T) xXx S'^-')^ 

\^\\w ■= ||'"||Li((0,T)xXx§'i-i) + 



, (2.12) 
(2.13) 



Li((0,T)xXxS'*-i) 



where ^ and f Vx are understood in the distributional sense. 

It is known [CI, C2] that the trace map 7_ (respectively 7+) from C"'^([0,T] x 
X X S-^-i) to C{X X S'^-i) X C{{0, T) X T_) (respectively C{X x S-^-^) x C((0, T) x r+)) 
defined by 

1-W = (V'(0,.)>^|{o,T)xr_) (respectively -f+{^) = (^^(T, .), V'|(o,T)xr+)) (2.14) 

extends to a continuous operator from W onto L^{XxE>'^'^ , t+(x, v)dxdv)xL^[[0, T) x 
r_ , min(T-t , t+ {x.v))dtd^{x., v) ) (respectively {X x S'^-^ r_ (x, v)dxdv) x ( (0, T) x 

r+, 

min(t, T_(a;, t'))(it(i^(x, f))). In addition 7-t admits a continuous lifting. Note that 
Li(X X S'^-i) is a subset of L'^{X x S'^-i, t+{x, v)dxdv). Note also that Li((0, T) x 
T_,dtd^) (respectively L^((0,T) x r+,(if(i(^)) is a subset of L^((0,T) x r_,min(r — 
t,T-^-{x.v))dtd^{x,v)) (respectively L^((0,T) x r+, min(t, t_(2:, t;))rft(i(^(x, v))). 

We now introduce the space 

W:={ueW\ 7-(n) G L^(X X S'^-i) X L^((0,T) x r_,rffrf^)}. (2.15) 

We recall the following trace results (see [CI, C2] in a more general setting). 
Lemma 2.1. The following equality is valid 

W = {ueW \ 7+(m) eL\Xx S'^-i) X L\{0, T) x r+, dtdC)}. (2.16) 

In addition the trace maps 

7± : ^ L\X X §^-^) X Li((0,T) x r±,dtdO 
are continuous, onto, and admit continuous liftings. (2-17) 



2.4 Solution to equation (12.21 ) 

For any r > 0, we identify the space i^^((0, r), L^(r-|-, d^)) with the space L}{{f), r) x 
T±,dtd^), and we extend any function G L^{{0,r), L^(T_,d^)) by on M\(0,r) 
(the extension is still denoted by 0). 

Let G Li((0,77),Li(r_,(iO)- We extend by outside (0,77). Then we 
consider the lifting G_(t)0 G of (0, 0) defined by 

G_ (t)0(x, t;) := e" .{x-sv,v)ds^_ _ ^_ ^) ^ ^ _ ^_ (^^ t;)^;, t;) , (2. 18) 
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for (t, X, f ) G (0, T) X X X S'* ^. Note that G^{.)(j) is a solution in the distributional 
sense of the equation + v'Vx)u + au = in (0,T) x X x E>'^~^ and 

||G'_(.)0||w < (l+l|o-||oo)||G'_(.)0||Li((O,T)xXx§d-i) < (l+||cr||oo)7'||0-||Li((O,r,)xr_,dtdO- 

(2.19) 

To prove the latter statements, one can use the change of variables given by Lemma 
4.1. From fl239D we obtain that the map i : L^{{0, rj) , L\T _ , d^)) W defined by 

z(0) = G_(.)0, (l>eL\{0,7j),L\T_,dO), (2.20) 

is continuous. 

The following result holds (see [DL, Theorem 3 p. 229]). 

Lemma 2.2. The equation i\2.2\\ admits a unique solution u in W which is 
given by 

u{t) = G-{t)(P+ [ U{t- s)A2G-{s)(t)ds. (2.21) 
Jo 

where U(t) is the strongly continuous semigroup in L^{X x S'^^^) introduced in sec- 
tion 2.2. 

From (12.201) . Lemma 2.2 and f l2.17p . we obtain the existence of the albedo 
operator. 

Lemma 2.3. The albedo operator A given by the formula 
A(p = M|(o,r)xr+, for G L\{0, 77), L^(r„, ^0) where u is given by (IXTTI) . (2.22) 
is well-defined and is a bounded operator from -^^-'^((0, 77), L^(r_, d^)) to L}{{f), T), L^(r+, d^)). 



3 Stability results for the inverse problem 
3.1 Recall of uniqueness results 

ChouUi-Stefanov [CSl] studied the uniqueness of the reconstruction of (a, A;) from 
the albedo operator by analyzing the distributional kernel of that operator. They 
considered the following problem 

Ou r 

— {t, x,v)+v'Vxu{t, x,v)+a{x,v)u{t, x,v) = / k{x,v' ,v)u{t, x,v')dv' , {t,x,v) eM^xXxV, 



dt 



(3.1) 



u\Rxr-it,x,v) = (j){t,x,v), 
M|i«o = 0, 

for (f) G Ljojjjp(M, L-'^(r„, d^)), where V is an open subset of M'', d >2. The albedo 
operator is defined as an operator from Ljj^j^p(M, L^(T_, d^)) to ^^^^^.(M, L^(r+, d^)). 
They proved, in particular, that the albedo operator uniquely determines the ab- 
sorption and scattering coefficient {a,k) provided that cr is a function of x and l^;! 
only. It is straightforward from the proof of this result (see [CSl, Theorem 5.1, 
Propositions 5.1 and 5.2]) that the following result holds. 
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Proposition 3.1. Assume that {a, k) are admissible and a{x, v) = f{x, \v\) for 
some real function f. Let T > rj > 0. Then the following statements are valid: 

iifT> diam(X) then the albedo operator 
A : L^{{0,r]), L^{r_,d^)) L^{{0,T),L^{r^,d^)) uniquely determines a, 

a ifT> 2diam(X) then the albedo operator 

A : L^{{0,rj),L\r_,d^)) L^{{0,T), L^{r+,d^)) uniquely determines {a,k), 

In this paper we analyze the stabihty of the reconstruction of (a, k) from the 
albedo operator. Our study is also based on the distributional kernel of the albedo 
operator. In a first stage, we do not assume that a{x,v) — f{x, \v\) for some real 
function /. 

3.2 Decomposition of the albedo operator 

Consider the distributional kernels 

a2(T,x,v,x',v') = / ^ '\-Io'^i^~P^,^)d^-Io-^'°~^^'^^'^i^-t^-P^'y)dP (3.3) 
Jo 

xk{x - sv,v',v)6:^-sv-t-{x-sv,v')v'{x')S{t - s - T-{x - sv,v'))ds, 

for a.e. (r, x, v, x' , v') G lRxr+ xr_ and where we have defined /gd-i fi{v')Sviv')dv' = 
Idx^<,(^')f^(^')(^l^(^') = M^'o) and J^6{T-s)fs{T)dT = /3(s) for {v,x'q,s) e 
S'^-i xdXxR and for (/i, /a, /s) G C{S^-^) x C{dX) x C(M). 

We consider the usual decomposition of the albedo operator as a sum of three 
terms: the ballistic part (whose distributional kernel is given by ai), the single scat- 
tering part (whose distributional kernel is given by 02) and the multiple scattering 
(whose distributional kernel is denoted by 0:3). Using [CSl, Theorem 5.1], we know 
that \u{x')v'\-^a3 e L^{r_, Ll^{R, L\r+,dC))). The following Proposition 3.2 im- 
proves on the latter statement provided that k e L°°{X x x S*^"-^). The result 
will be used in the proof of Theorem 3.1. 

Proposition 3.2. Assume d >2 and {a, k) admissible. Assume that k G L°°{X x 

A(0)(t, x, := / [ai + a2 + OL^){t — t\x,v,x' ,v')(t){t' ,x' ,v')dt' diJi{x')dv' , (3.4) 
i(o,r,)xr_ 

for a.e. {t,x,v) G (0,T) x r_|_ and for any continuous and compactly supported 
function on {0,r)) x r_, where 

\u{x')v'\-'a3 G L-(r_,i7((-r7,r),i7'(r+,c/0)), for any 1 < p < (3.5) 
Proposition 3.2 is proved in section 4. 
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3.3 First stability result 

Now we assume that X is a bounded open convex subset of R*^, d > 2, with 
boundary and that 

the function < a is continuous and bounded on X x S*^^^, 



the function < /c is continuous and bounded on X x S'^ ^ x S'^ ^. 



(3.6) 



Let {a, k) be a pair of absorption and scattering coefficients that also satisfy 
( KB . Let i be the albedo operator from L^{{0,r]), L\T_, d^) ) to ( (0, T) , {T+ , d^) ) 
related to (a, A;). 

For {x,v,s,w) e r± X M X S'^"^ < s < t^{x,v), let E^{x,v,s,w) > be 
defined by 

(PS pT^{x^SV,w) 
— I a{x^pv,v)ds— / a{x ^ sv ^ pw,w)dp] . 

(3.7) 

Replacing cr by cr in f l3.7l) we define Ezf{x,v,s,w) similarly for {x,v,s,w) G F-i- x 
M X S*^-!, < s < T^ix,v). 

Let {x'q, v'q) e r_. For £i > and 62 > 0, let fe, E C\T.) and G C°°(M) be 
such that 

Ai > 0, supp/e, C {(x', v') eT_ \ \x-x'q\ + \v' -Vq\< £1}, (3.8) 
fe,{x',v')d^{x',v') = 1, 



r_ 

/•+00 

5'£2 > supp5fj2 ^ (0,min(?7,£2)), / Qe^dt = 1. (3.9) 

Consider the function 0ei,£2 G C^(M x F_) defined by 

0e,,e2(^',a:',v') = ge2it')fe,ix',v'), (3.10) 

for t' G (0,+oo) and ix',v') G F_. Note that supp^e,,^^ C (0,?]) x F_ (see (K9h). 
From (13.81) and (13.91) it follows that \h'{x')v'\(psi,e2 is a smooth approximation of the 
delta function on M x r_ at (0, Xq, v'q) as Si 0+ and £2 — ^ 0+. 

Let ip be any compactly supported continuous function on (0, T) x F+ such that 
1 1 1 1 00 < 1- First we note that upon using the estimate ||'0||oo < 1 and the equality 
/(Or?)xr 4'ei,e2iti ^^''^)dtdC,{x,v) = 1 wc obtain that 

< II - l|Li({0,T),Li(r+,dO) 

< P-^lk,r, (3.11) 



^p{t, x,v) [A — A) 4>£i,£2it, X, v)dtd^{x, v) 

(o,T)xr+ 



where ||.||,,,t := || •IU(Li((o,»?),Li(r_,dO),Li{(o,T),Li(r+,dO))- 

In addition, it follows from (I3.4l) - (l3.5p that for any compactly supported con- 
tinuous function -0 on (0, T) x r+, we have 



ij{t,X,v) [A~ A) (j)£^^£^{t,X,v)dtd^{x,v) = Ji(^/^, £i, £2)+/2(^/', ^l, f^2)+-^3(V', £i, ^2), 

{o,T)xr+ ^ ^ 

(3.12) 
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where 

x?/j(t, X, f )0ei,e2('^ ~ v),x — r_(x, v)v, v)dC,{x, v)dt (3.13) 

/2(^,£i,£2) = / / '^{t,x,v) / / (E_{x,v,s,v')k{x - sv,v',v) 

Jo Jr+ Jsd-i Jo ^ 

—E-{x, V, s, v')k{x — sv, v', v)j 'ipei,e2{'t — s — r_(z — sv, v'), 
X — sf — r_(x — sv, v')v' , v')dsdv'd^{x, v)dt (3-14) 



|/3(^,£i,e2)| < m,x,v)\Pdax,v)dtj , (3.15) 

where C is a constant that does not depend on Si, 62 (for fl3.15p we also used Holder 
inequality and the equality \\(t)e^,e2\\L^{{o,r,),L^(T.,dO) = 

Using fl3.8p - fl3.10p . fl3.6p . we obtain the following preparatory Lemma 3.1. 

Lemma 3.1. Assume that X is convex and that {a, k) and (a, k) both satisfy (13. 6p . 
Then the following statements are valid: 

i. ifT> diam(X) then 

lim Mm hiip, 61,82) = ipiT+{xQ,VQ),XQ + T+{xQ,VQ)v'f^,VQ) (3.16) 

£2— >0+ £1^0+ 



X 

for any compactly supported and continuous function if) on (0,T) x r+; 

ii. ifT> 2diam(X) then 

lim lim 12(^,61,62) = I^W + lUi^), (3.17) 
£2— >-o+ £1^0+ 

for any compactly supported and continuous function ip on (0,T) x r_|_, where 

Ilii') = / Ipis + T+{Xq + SVq,v),Xq + SVq + T+{Xq + sVq,v)v,v) 

{k — k){x'i^ + sfg, v'q, v)E^{x'q, v'q, s, v)dsdv, (3.18) 
J2W = / ip{s + t+{xq + sv'q, v),x'q + sv'q + T+{Xq + sVq,v)v,v) 

J§d-1 Jo 

k{x'Q + sv'q, v'q, v){E^ — E^){x'q, v'q, s, v)dsdv, (3.19) 
where E+ and E^ are defined by (13. 7p . 

Lemma 3.1 is proved in section 5. 

Taking account of Lemma 3.1 and (13. lip , and choosing an appropriate sequence 
of functions "?/'" , we obtain the main result of this paper: 

Theorem 3.1. Let T > r] > 0. Assume that d > 2 and X is convex and {a, k) and 
(a, A;) both satisfy condition (13. 6p . Then the following statements are valid: 
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i. ifT> diam(X), then 



exp ( — / cr(xQ + sVq, v'Q)ds 



— exp 



a"(xQ + sVq, v'Q)ds 



(3.20) 



a. ifT> 2diam(X), then 



S'^-i Jo 



k-k 



{x'q + sv'q, v'q, v)E+{xq, v'q, s, v)dsdv 



<T+{Xq,Vq) sup Upy^Q 

se(o,T+(x^,i>^,)) 



q,uq) sup 

3S(0,T+(a:{,,u^)) 
t,gSd-l 



{x'q,v'q,s,v) 



+ U-A\\ 



(3.21) 



where ||.||^,t := ||-||£(Li((o,r,),Li(r_,de)),Li((o,r),Li(r+,dO)) "'^^ where E+ and E+ are de- 
fined by (1X71) . 

The proof of Theorem 3.1 is given in section 5. 

Remark 3.1. One can prove that estimate (13.201) still holds 
provided that T > diam(X) and k G L°°{X x S"'^^ x S'^^^) and without assuming 
(13.61) nor that X is convex. 



3.4 Second stability result 

We now impose that the absorption coefficient a does not depend on the velocity 
variable, i.e. a{x,v) = a{x), x & X. Then let 

M := {{a{x),k{x,v',v)) e L°°{X) X L^{X xS'^-^ x^''-^) \ {a,k) satis&es 

and a E H^+^iX), lk||^.+,^^^ < M, < M}, (3.22) 

for some f > and M > 0. Using Theorem 3.1 for any {xq, v'q) G r_ we obtain the 
following Theorem 3.2. 

Theorem 3.2. Assume that d > 2 and X is convex. Let T > r] > 0. For any 

(a, k) E M. and (a, k) E M. the following stability estimates are valid: 

i. ifT> diam(X) then 

\\a-a\\Hsix)<C,\\A-A\\lT, (3.23) 
where -1 < s < | + f, k = '^j^^, and Ci = Ci{X, M, s, r); 
a. ifT> 2diam(X) then 

k{xQ + s'vq, v'q, v) — A;(xg + s'v'q, v'q, v) ds'dv 
< C,\\A - Ar,^^ [l + \\A - A\\l^) , (3.24) 

for {x'q, v'q) E r_, and where k = j^^^^, < r < f, and C2 = C2{X, M, r, r); 
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in. in addition, if T > 2diam(X) then 



\\k - A;|Ui(xx§<*-ixs<^-i) < CM - A\\l^ (I + \\A - iHj^/j , (3.25) 

where k = -^^fj^, < r < r, and C3 = Cs{X, M, r, f). 
Theorem 3.2 is proved in section 5. 

Remark 3.2. Stability estimates similar to (13.231) were given by Cipolatti- 
Motta-Roberty [CMR, Theorem 1.1]. They proved (I3.23P for s = — | under the as- 
sumptions k,ke L°°{X, L^(§'^~-^xS'^" -*-)), max(||a||oo, ||5"||oo) < M (and max(||crp||oo, ||5'p||oo) < 
00). They also proved (13.231) for < s < f under the assumptions k, k E 

L-^iX,L\S^-^ X S'^-i)), a,ae Hi+^\X) and max(||a||^.+„ M^^^,) < M. 

4 Proof of Proposition 3.2 

Before giving the proof of Proposition 3.2, we need Lemmas 4.1, 4.2, 4.3. 
Lemma 4.1. For f E L^{X x E>'^~^) we have 

f{x,v)dxdv= I I f{x^sv,v)dsd^{x,v). (4.1) 



xxs-^-i Jr± Jo 



For the proof of Lemma 4.1, see [CS2, Lemma 2.1]. 

Let m > 1. For U a subset of R™', we denote by xu the function from R'" to 
defined by 

1 if X eU. 



Xu{x) - |q otherwise. ^^'^^ 

Lemma 4.2. Let T > and let 1 < p < Consider the nonnegative measurable 
function /9 : X — R defined by 

T+diam{X) _ (a; _ a;')^^p('i~3) 



y{d~ 



(3(x')= / / ^ ^ '--^T^:r^dsd^(x,v) (4.3) 

for a.e. x' E X . Then 

P e L°°(X). (4.4) 



Proof of Lemma 4-2. We first consider the case d = 2. 
We have 

/3{x') < [ r [ 7 , ^ „ dvdsd^ix) 



< 



dxJ\x-x'\Js^ [S- [X-X')V)P 

2n r+00 _l I 

Tl~7 \ 7\ —^dujdsdfi{x) 

dx Jo J\x-x'\ p-lds{s-\x- x'\ cosuj)P 

1 /"^^ 1 

dx / — ^ — rduj. 



p — 1 Jqx \x — x'\P~^ Jq (1 — cosa;)P~^ 
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for E X. Hence using the estimate p < 1 + |, we obtain 

Now assume d = 3. Using fl4.3l) . spherical coordinates and performing the 
change of variables "s" = |a; — we obtain 

JS2 J\x-x'\ \SV — [X — X )\ f 

= 2n — -dii(x) I — -. / ^ — -duds 

Jqx \x - x'|2p-i ' 2s{p - 1) duj (s2 + 1 - 2ssmuj)P-^ 

f 1 f°° I 
< 271 — -diiix) / — tt; — -ds. 

Therefore using the estimate l<p<l + |, we obtain 

1 , . . r 1 



||/?|U.(^)<2vrsupy^^^-^rfM-)y^ 1)2.-. ^^ (4-6) 

Finally assume (i > 4. Note that |s— = — (x— | < \sv-~{x—x')\ 
for s G R and x, x' G M'^. Using in particular the latter estimate and (14.31) . we obtain 

{p{d - 1) - 1) Ix £-1 (11^ -x'\v-{x- x')\yi^-^)-^'^'''^^^''^ 



1 . , . /" 2 COS w'^ ^ 



for G X. Hence using the estimate p < 1 + we obtain 

1 r 1 



L-(x) < Vol(§ )--— -— — -iciiTprsup / 7— -^(^ZTpT^^(^) 
{p{d- I) -1)2 — 2 — ^exJaxF-^^r^ ' 

2 coso;'' ^ , ,^ _ 

. < oo- (4-7) 

-f (1 — sinw) 2 

□ 

Finally, we need the following Lemma 4.3. 

Lemma 4.3. Consider the nonnegative measurable function 7 : (0, T) x X x S'^^^ x 
X X §"'-1 ^ M de/ine(/ &?/ 

7(t, X, V, X', V') = 2'^'\^o^t)i\x - X'\) [e"^o' a(x-sv,v)ds-J^~^^ a(x-s^v-pv^,v^)dpQ^^^ ^ _ ^^^^ 

X 6{x - siv,x')k{x - siv,vi,v)k{x',v',vi)] _ t2_(,_,/)2 ^ — — — ^ , , (4.8) 

''^-2(t~(x-xl)v) \tV — X — X'\^'^ * 



where 6 is defined by (I2.10p . Then 



Ui{t — si)A2Ui{si)A2fdsi\{x,v)= / ^{t,x,v,x ,v)f{x ,v')dx'dv 

(4.9) 

fort G (0,r) and for a.e. {x,v) G X x S'^-i and for f e L\X x E,'^-^). 
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Proof of Lemma 4.3. Let t G (0, T) and let f e L^{X x S'^-^). From ([23]) and (1231) . 
it follows that 



'u,{t-s,)A2U,{s,)A2fds, = 11 e 




/q (j(x—pv^v)d'p—^Q^ a{x~{t—si)v~pvi,vi)dp 



'0 JS<*-ix§d-i 

X — (t — Si)!), Vi, t;)/i;(a; — (t — Si)v — SiVi, v' , Vi) 
x6{x — — si)v, x)9{x — (t — si)v — siVi, x — (t — si)!;) 
xf{x — {t — si)v — siVi, v')dv' dvidsi. (4-10) 

Performing the change of variables "si = t — si" and then performing the change 
of variables 'V = x — {t — si)vi — siv" (2*^"^ ^''^i^l^^/I^^pdL dx' = dvidsi), we obtain 

iD. " " □ 

Proof of Proposition 3.2. Let G L^((0, L^(r_, d,^)). Let u be the solution of 
(12.21) . Using twice Duhamel's formula (12. lip and using (I2.2ip we obtain 

u{t) = Ri{t) + R2{t) + Rl{t) + Rl{t), (4.11) 

for t G (0, T) where 

i?i(t) = G4t)<j), (4.12) 

R2{t) = [ Ui{t-t')A2G-{t')4>dt', (4.13) 

Rl{t) =11 Ui{t-t' - si)A2Ui{si)A2G4t')(f)dsidt', (4.14) 
Jo Jo 

Rlit) = / U^{t-t' -S2-Si)A2U^{s{)A2 (4.15) 

Jo Jo Jo 

U{s2)A2G.{t')(f)dsids2dt'. 
From (14121) and (EH]), it follows that 

^i|(o,T)xr+(t, X, v) = e-Zo^"'""' - r_(x, t;), x - t_(x, v)v, v) 

= / ai(t-t',x,t;,x',f')0(t',x',t;')c?t'(i^(x',f'), (4.16) 
i(o,»7)xr_ 

where ai is defined by (13. 2p . 

From ([23]), fl2J8l) and (jUS]), it follows that 

R2{t,X,v) = / e{x-t'v,x) k{x-t'v,v',v)e-^0 -^^-P^'''^'iP-fo~ ^ a(x-t'v-pv',v')dp 

Jo i§d-l 

x^(t — t' — r-{x — t'v, v'),x — t'v — r_(x — t'v, v')v' , v')dv'dt' . 

Hence 

-R2|(o,r)xr_ 3;, t;) = / a2{t — t' , x,v,x' ,v')(p{t' , x' ,v')dt'd^{x' ,v'), (4.17) 

i(o,J7)xr_ 

where 0^2 is defined by (13. 3p . 
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From (gUD and (gS]) (with "f = t-t'), it follows that 
Rlit) = / X(o,+oc){t-t')j{t-t',x,v,x',v')G^{t')(f){x',v')dt'dx'dv', 

i(0,+oo)xXxS''-i 

(4.18) 

where 7 is defined by (14.81) . Using (I4.18P and (14.81) we obtain 
RlifQT)xT {^1^1'^) = I ds^{t — t' , x,v, x' ,v'){G_{t')(l)){x' ,v')dt'dx'dv' ^ 



where 



a3^{T,x,v,x',v') = X{Q,+oo){r - \x - x'\) 



(4.19) 



■ Jq^ cr{x—pv,v)dp— Jq a-{x — siv—pvi,vi)dp 



X k{x — Siv, fi, v)k{x', v', Vi)9{x, x — Siv)9{x — Siv, x')] _ ^2_|^_^/[2 

"''-^ 2(t-Ix-x')v) 



x — x — s-i V 



x2 



d-2 



[t — {x — x')vY ^ 

\x — X' — Tf P'^"^ 



(4.20) 



for a.e. (r, x, v, x', v') eRxT+ x X x S"'"^ 
From fHTTHD and it follows that 



Rlit) 




t rt-t' p 
Jo ixxSt*-! 

where 7 is defined by (14.81) . Hence 
-^3|(o,r)xr+('^' ^' = / / / 

Jo Jo JXxSd-i 

for a.e. {t,x,v) G (0,T) x r+ where 
^{r,x,v,x',v') = 2'^~2x(o,r)(|a; - x'\) 



^{t—t'—S2, X, V, x' , v'){U{s2)A2G-{t')(/)){x' , v')dx'dv'ds2dt' , 

(4.21) 



^{t-t'-S2, X, V, x', v'){U{s2)A2G^{t')(p){x', v')dx'dv'ds2dt', 

(4.22) 



- cr(x—sv,v)ds—jQ a(x—siv—pvi,vi 



X,X — SiV) 



X 6{x — SiV, x')k{x — SiV, Vi, v)k{x' , v' , Vi)] _ ^2_(^_^/)2 



(r — (x — x')v) 



f\„,\d-3 



' 2{r-{x-x')v) \TV — X — X' 



l\2d-i 



(4.23) 



for a.e. (r, x, v , x', f ') G (0,T) x r+ x X x S'^-i. 

Let i) G L°°((0, T) x T+). Assume that k e L°°{X x S^^-^ x S'^-^). From (ICTll 
it follows that 



r-T 

/o 4 




?/^(t, X, t^)i?3|(o,T)xr+ ^' v)dtdi{x, v] 



(G_(O0)(x',t;') / / a~3'(t-t' ')4'{t, X, v)d^ {x, v)dtdt'dx'dv' 

(o,r)xXxS''-i Jo Jr+ 



< ||G'_(.)0||ii((o,T)xXx§d-i) 




«3 {t — t',x,v, x', v')il){t, x,v)d^ {x, v)dtdt' 



JV+ 



L-(IR,,xX,,x§f,-') 



(4.24) 



13 



From Lemma 4.1 and fl2.18p . it follows that 

\\G-{t')cf)\\LHXxSd-^) < ll0l|Li((O,r;),Li(r-,rf?))' ^' ^ i^^^). 

From fl4.20p . Holder's inequality and Lemma 4.2, it follows that 



(4.25) 




tts (t — t' , X, V, x' , v')'i/j{t, X, v)dC,{x, v)dt 



JT^ 



< 




JV+ 



< 2 



d-2\\ T ||2 



X(o,+oo)(^ -t' -\x- x'\){t -t' - {x- x')v) 
2^-'^\x-x' - {t - t')v\'^'i~^ 

\ij{t,x,v)\P'd^{x,v)dt] (3{x')v 



-\ip\it, X, v)dC,{x, v)dt 




JT+ 




<2'^-1|fc||^||/?|II.(^) / / m,x,v)\^ di{x,v)dt 



Jv^ 



(4.26) 



for a.e. {t',x',v') G (0,T) x X x S'^-i. 
Using flT2D -f l4^ we obtain 



r-T 




i){t, X, v)-R3|(o,T)xr+ 2;, v)dtd^{x, v 
<2''-^T||fc||L||/?||f^(^)||0|Ui((o,,),Li(r_,dO) 



'0 Jr+ 

•,d-2rp\\ T ||2 



(4.27) 



i-T 




\'4j{t,x,v)fd^{x,v)dt 



Jr. 



In addition from fl4.22p . Holder inequality, fl4.23p and Lemma 4.2 it follows that 



i-T 




Jr+ 



ip{t, X, t^)^3|(o,T)xr+ 2;, ?;)rf^(x, f )dt 




t rt-t' 



Ipit, X, v) 



Jr. 




JO 



7(t — t' — S2, x, V, a;', v') 




Jr+ 



X(t'+s2,T)(^)^(i,2;,'y) 



X (t/(s2) AaG, (t')0) (a;', v')dx'dv'ds2dt'd^{x, v)dt\ 

If [ iUis2)A2G4t')<p)ix',v 
X 7(t — t' — S2, X, V, x', v')dtd^{x, v)dx'dv'ds2dt' 

< r r ' [ \u{s2)A2G4m{x',v') 

Jo Jo JxxS'i-'^ 

i / X{t'+s2,T){t)\"y\^i't — t' — S2, x,v, x' ,v')dtd^{x,v) ^ ds2dt'dx'dv' 

JV+ 



\ip{t, X, f)|^ dtd^{x, v] 



(o,T)xr+ 



<2'^-^||fc||Lll/?ll^ 




'0 -'0 JXxnd-'^ 

\i{j{t, X, v)\^ dtd^{x, v) 



\U{s2)A2G-{t')(t)\ {x',v')dx'dv'ds2dt 



{o,T)xr+ 



(4.28) 
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Moreover using (14.251) . the equality 11^2 1 1 = ||c'"p||oo and the estimate ||t/(s2)|| < 



e^all'^pll- for S2 > (see Trotter's formula [T] [/(sa) = s - lim„^oo )e 
where Ui and A2 are defined by (12. 9p and (12. 6p respectively), we obtain 

\U{s2)A2G4t')(f)\ {x',v')dxdv'ds2dt 




'0 ^0 



< 




^0 



dS2 



Li((0,r))xr_,dtdO 



< T{e 



T\Wp\\ 



Li((o,7?),Li(r_,dO)- 



Combining (KT7\\ - K^ . we finally obtain 




Jr^ 



<2'^-lA;||L||/3||I^(^)Te^ll'^^ll-||0|Ui((o,,),ii(r_,.o) 

1 

|'0(t, X, f)|^ (i^(x, v)dt 








Proposition 3.2 follows from fl416|) . fliTTD and flCTj) . 



(4.29) 



(4.30) 
□ 



5 Proof of Lemma 3.1, Theorems 3.1, 3.2 

Proof of Lemma 3.1. First note that using twice Lemma 4.1 we obtain 

I I f{x—wv,v)dwd^{x,v)= / f{x,v)dxdv= / / f{x'+sv',v')dsd^{x',v') 
Jr+ Jo ixx§d-i ir_ Jo 

(5.1) 

for / G Li(X X S'^-i). 

We first prove (13.161) . Let T > diam(X). We have, in particular, T > t_{x,v) 
for any {x,v) e X x S^-^. From (1^71^ and fIXTUD it follows that 

/i(?/',£i,£2) = / $£2(x, v)/£,(a; - r_(x,t;)t;,t;)rf^(x,i;), (5.2) 



+ 



where is the continuous function on r_|_ given by 

^,^{x,v) = [ ij{t,x,v)g,,{t-T_{x,v))dt (e- ^o'^"^^ _ lo-^'''^^ H--sv,v)ds 



Jt—{x,v) 

(5.3) 

for {x,v) G r+ (we used also supp^f^j C (0,+oo)). The continuity of follows 
from the assumptions (13.61) . the continuity of tf) and and the continuity of r_ 
on r+ (X is convex with boundary). From (15.11) ("/(x,f) = ^^ $^2(3; + 
T^{x,v)v,v)fe2{,x — T^{x,v)v,vy') and (13.81) . we obtain 

^e2{.^^'").feAx - r-{x,v)v,v)d^{x,v) = I $£2(3;' + T+{x',v')v',v')fe^{x',v')dC{x',v') 

r+ Jr- 

xg,^{t-T+{x'Q,v'Q))dt. (5.4) 
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The limit (l316|) follows from (JED, dEl, and the continuity of and (KB . 

We prove (13.171) . Let T > 2diam(X). We have, in particular, T > s + t^{x' + 
sv',v) for s G (0, r+(x', f')) and (a:',^') G r_. From (13.141) , it follows that 

P /•T-{x,v) pT 

I2{^,£i,e2)= / / / i'{t,x,v) (5.5) 



JT+ Jo Jo 

X / {k{x — wv,v' ,v)E^{x,v,w,v') — k{x — wv,v\v)E_{x,v,w,v')) 
^9e2{t — w — T_{x — WV, v'))fsi{x — WV — r_ (x — wv, v')v' , v')dv'dtdwd^{x, v). 
Using (15.51) and (15. ip . we obtain 



r_ 



/2(^,£l,£2)= / ^e2ix',v')f,,ix',v')dax',v'), (5.6) 

where 



p rr+{x',v') pT 

'^e2ix',v')= / / / 'ilj{t,x' + sv' + T^{x' + sv',v)v,v) 

JS'*-l Jo J s+T+(x'+sv' ,v) 



X yk{x' + ', f ', v)E^{x' , v', s, v) — k{x' + sv', v', v)E^{x , v' , s, v) 
y<ge2{t — s — T^{x' + sv' ,v))dtdsdv, (5.7) 

for {x',v') G r_. From (13. 6p . (13. 7p and the continuity of ip and gs2, it follows that 
is continuous on r_. From (13.81) and (15. 7p it follows that 

'^e2ix', v')fe^ {x', v')di{x', v') 
-r+ix'o^v'^) ( pT 



%l){t, x'q + sv'q + T+(xo + sv'q, v)v, v) 

ei^O+ Jgd-i Jq yj s+T+{x'g+sv'g,v) 

(E+ix'Q, v'q, s, v)k{x'Q + sv'q, v'q, v) - E+{x'q, v'q, s, v)k{x'Q + sv'q, v'q, v)) (5.8) 
X ^(^ "~ -5 — '^+{^0 + ■^''^0' 'L'))(^^) dsdv. 

The limit (IHTTI) follows from ([52]), ([53]), ([S3]), ([SS]), the continuity of ip and the 
Lebesgue dominated convergence theorem. □ 

Proof of Theorem 3.1. We first prove Let T > diam(X). Let £3 > and let 

ijj^^ be a continuous and compactly supported function on (0, T) x r+ that satisfies 

< V^ea < 1 and suppz/'gg C {{t, x, v) G (0, T) xT^ \ \v - v'q\ < £3}, (5.9) 
'0e3(t,x,f) = 1 for {t,x,v) G (0,T) x r+ such that 

1 ' I ^ ^3 1+ / / '\l ^ ~ ''"+(^0' ''^0) 



2" -^v-u^^u/i- 2 
From fl316D and ([ES]), it follows that 



/i(^.3,^i,£2) = e-/o^^*^°""'-(^o-H.'^i)'^^ _ (5.10) 
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Using (13.171) . (15.91) and the estimate a > 0, we obtain 



lim lim /2(V^e3,£i,£2) 

£2^0+ £1^0+ 



<diam(X)(||A;|U + ||^||oo) / ^ dv. 



Hence 



lim lim lim /2 (■?/'£.., £i, £2) = 0. (5-11) 

£3^0+ e2-»0+ £1-^0+ 



From (13151) and (ESD, it follows that 

'0 JdX 



/3('0e3,£i,£2)| < Ci I I I dvdti{x)dt 

^Jo J dx J .-"'^''f;^^ 



for £j > 0, i = 1 . . . 3. Therefore 



lim limsup limsup J3('0e3, £1, £2) = 0. (5-12) 

£3^0+ £2^0+ £1^0+ 

In addition, from (l3TTil - ( 13121 ) it follows that 

|/l(?/^e3,£l,£2)| < \\A-A\\^^T+ |/2(^£3,£l,£2) + /3(V^£3,^1,£2)|, (5.13) 

for £j > 0, i = 1 . . . 3. 

Combining (15101) - (l57[3l) we obtain (Km . 

Now we_ prove (jMH). Let T > 2diam(X). Let U := {it',v) G {0,t+{xq,v'q)) x 
§"'-1 I {k - k){x'f) + t'v'f^,v'f^,v) > 0}. From 1^ it follows that U is an open subset 
of M X S'^^^. Let (Km) be a sequence of compact sets such that UmeN -^"i = U 
and Krr, C Km+i for m G N. For m G N let Xm e C°°(M x §'^'\R) such that 
XKrn < Xm < Xu (where xk^ and xu are defined in (H2])), and let 

Pm = 2xn.-l. (5.14) 

Thus we obtain 

lim {k-k){x'Q + t'vQ,v'Q,v)prn{t',v) = \k-k\{xQ + t'v'Q,VQ,v), (5.15) 

m— >+oo 

for V G S'^-i and t' G (0, r+(x;„ i;^)). 
Consider 

V5:={(t,x,t;) G (0,r) xr+ I \v - {vv',)v',\ > 6, ^-<t<T-^-}, (5.16) 

V5,i := {(t, a;, z;) G (0, T) x r+ | - (OT(,)t;(,| > 5 + y , 5 < t < T - 5}, (5.17) 

for < 5 < and / G N, ; > 2. For < 5 < 2^^^^ and / G N, / > 2, let 

be a continuous and compactly supported function on (0, T) x r+ such that 

XVs, < X5,i < XVs (5.18) 
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(where Xvsi XVs defined in (14.21) ). Finally, for < 5 < and m,l G N, 

/ > 2, let ■05,m,/,£3 be the continuous compactly supported function on (0,T) x r+ 
defined by 

tljS,ra,l,e3{t,X,v) := X5,lii^^^^)iCe-sit- S- s')pmis,v)) , (5.19) 

where Ces e C^(M), Ce^is") = 1 for s" e [-63,63], < Cea < 1 and Ce,{s") = for 
\s"\ > 26s. 

From (I3.16p . (I5.19P and the equality x<5,/(^) ^^o+'^'+l^O' ^o)^0' "^0) = for t G (0, T) 
(see (l5l6l) -( l5:T8ll ). it follows that 

lim lim Ji(V'5,m,/,£3, £1, £2) = (5.20) 
for < 5 < HHBiiZ), m,l>2,6s> 0. 
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From (1318I) -( 1339I) and (ICTl) . it follows that 
Jsd-i Jo 

xpm{s, v){k — k){xQ + sv'q, v'q, v)E+{xq, Vq, s, v)dsdv. (5.21) 

l2ii'S,m,l,e3) ■= / Cesi0)X5,li^ + T'+i^0 + Sv'q,v),Xq + SVq + T+{x'q + SVq,v)v,v) 

xpmis, v)k{xQ + sv'o, Vq, v) (^E+ - (xq, Vq, s, v)dsdv, (5.22) 

for < 5 < m, / > 2, £3 > 0. 

Note that using f l5.16l) -( l5T8l) we obtain 

lim XVs I (t, X, v) = xvs (t, X, v) , (5.23) 

t— >oo 

for {t, X, v) G (0, T) X r+ and < 5 < HliEOZi . 

From equality CeslO) = 1, (15.211) . (15.231) and the Lebesgue dominated conver- 
gence theorem, it follows that 



,lim /2('05,m,/,£3) = / / xvsis + t+{xq + sVq,v),Xq + sVq + t+{xq + sv'q,v)v,v) 

Jq 

xpm{s,v){k - k){xQ + svQ,v'Q,v)E+{xQ,VQ,s,v)dsdv, 



for < 5 < m G N, / > 2, £3 > 0. Therefore, using (15.151) and the Lebesgue 

dominated convergence theorem, we obtain 

lim lim /2^(V^5,„,i,£3) = / / xvsis + ^+{^0 + sv'q,v), x'q + sv'q + t+{xq + sVq,v)v,v) 



x_E'_i_(xo, Vq, s, v)\k — k\{x'Q + sVq, v'q, v)dsdv, 
for < 5 < EHB^LiZl^ £3 > 0. Using this latter equality and (I5.16p . we obtain 

lim lim lim lim l2iips,m,i,£s) = / / \k—k\{xQ+svQ,VQ,v)E+{xQ,VQ,s,v)dsdv. 

<5^0+ £3^0+ m^+oo i-»+oo ' '' J§d-i Jo 

(5.24) 
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From equality Cea (0) = 1 and (15.221) it follows that 
|/2(V'5,m,i,e3)l < diam(X) sup ap{xQ+sVo,VQ) sup {E+ - E+){xq,Vo, s,v) 

(5.25) 

for < 5 < m G N, / > 2, £3 > 0. 

Note that using (15.191) and the estimate < pm < 1 for all m, we obtain 



\lpS,m,l,e3{t, X, V) |P d^(x, v)dt 

< r I X.X^,^,-r'C3(t-^-^0^_u..-..,., dax,v)dt, (5.26) 



Jr 

T 



>''r+ l-(„;^)2 



for < 5 < EHH^, m G N, / > 2, £3 > 0. Therefore using the definition of we 
obtain 

\ij8,m,l,S3{t,X,v)fdC{x,v)dt< / Ix[~2ez,2ei\{t-S-S ) (i^(x, 

r I " 



^a}^a\>^ ^ — T — '1 / s'> 



(5.27) 

for < 5 < m G N, ^ > 2, £3 > 0. Using (15:261) - (15:271) and the Lebesgue 

dominated convergence theorem, we obtain 

lim limsuplimsup / / \ips^rn,i,e3{t,x^v)\^ d^{x,v)dt = 
£3^0+ „i-*+oo Jo Jr+ 

for < 5 < EHBLLiU, Using this latter equality and (13.151) . we obtain 

lim lim lim sup lim sup lim sup limsup |/3('05 m,/,£3 5 ^i, ^2)! = 0. (5.28) 

(5^0+ £3^0+ 1^+00 £2^0+ £1^0+ 

In addition, from (I3.1ip - (l3.12p . it follows that 

)|, (5.29) 

for < 5 < m EN, I >2, Ei > 0, i = 1...3. Combining (15:291) . (15:241) . 

(K^ and (IS2HD we obtain (^M)- □ 

Proof of Theorem 3.2. The method used to prove ( I3.23P is the same as in [W] and 
[BJl]. For the reader's convenience, we adapt the proof given in [BJl] with minor 
modification. 

Let (cr. A;), (a, k) E M.. We extend a and a outside X by 0. Let f = a — a and 
consider Pf the X-ray transform of / = a— a defined by Pf{x, Lp) := f{tip+x)dt 
for (x,v?) G TS'^-^ := {{z,v) G M"^ x S'^'^ \ vz = 0}. 

From fix G hI+^'{X), it follows that 

^-i(^)<I>i(rf,X)||P/||„ (5.30) 
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where 



11^/11*: = 



\Pf{x, ip)\'^dxd(p 



and Di{d,X) is a real constant which does not depend on / and II;^ := {x G 
R"^ \ xifi = 0} for if G S'^-^ Note that Pf{x,ip) = for (x, (/?) G TS'^"^ and 
|x| > sup^gj(^ \z\. Therefore using also (15.301) we obtain 

||/||^-i(^)<I^2(rf,X)||P/|U«>(^g.-l), 

where D2{d, X) is a real constant which does not depend on a, a. 
We also use the following interpolation inequality: 



(5.31) 



d+2f 



H^{X) 



< 



'^^tr 11/11 rr, 



for -i < s < f + r. As (cr, A;) G A^, it follows that 



\a\U < Ds{d,f)\\a\\^,^, < D3id,f)M. 



Therefore, 



<j{xq + sfo)(is < dia.m{X)D3{d, r)M, 



(5.32) 

(5.33) 
(5.34) 

\P{a-a){xQ,Vo)\, 
(5.35) 

for {x'q, Vq) G r_ (we used the equality e*^ — e*^ = €^^(^2 — ^i) for ti < t2 E M. and for 
some c G [^1,^2], which depends on ti and t2)- 

Combining f03|) . flOTD and flX^ . we obtain 



for (xQjfp) G r_. From f l5.34p . it follows that 



> e 



-diam(X)Z)3((i,f)Af 



-diam(X)D3 ((i,f)M 



(5.36) 



Combining fOUD and ^M), we obtain flX^ . 

We now prove (13.241) . Using (13. 7p and (15.331) . we obtain that 



k{xQ + sv'f), Vq, v) — k{xQ + sv'q, Vq, v) E{xq, Vq, s, v)dsdv 



S'^-i Jo 

> g-2diam(X)D3Kf)M / / \k - k\{x'Q + Sv'q, v'q, v)\dsdv , (5.37) 

Jsd-^ Jo 

for any {x'q,v'q) G r_. 

As {d-,k) E A4 we have ||5"p||oo < M. Using the latter estimate and (13.71) . we 
obtain 



sup \E - E\{x'q, v'q, s, v) < Me=^'^'^"'(^)^3('^'")*^ 



sup crp{xQ + sv'q, v'q 

seiO,r+ix'o,v'o)) (4.^(,)er_ 

X sup l\a — a\{x — pv,v)dp + / \a — a\{x'Q + sv'q + pv,v)dp 

j:'f,,v'f,)eT_ Jo Jo 

3e(0,T+(3:J,,t,^)) 
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< 2diam(X)Me2'^'"'"(^)^^('^'")^^||(T - a\\^, (5.38) 

for any {xq,v'q) E r_. (We also used |e" — e"| < e™'^^^'"'''"'-'|M — u\ where u = 
— a{xQ+pVg, VQ)dp — f^'^^^"'^^^"'^^ a{x'Q + svQ+pv, v)dp and u denotes the real num- 
ber obtained by replacing o" by 5" on the right-hand side of the latter equality that de- 
fines u ; using (15.331) (for a and for a) we obtain max(|tt|, \u\) < 2dia.m{X) D^^d, f)M.) 
Note that ||cr — 5"||oo < D3{d,r)\\a — 5'\\^d_^_^ for < r < f (see (15.331) ). Therefore, 
combining (15371) . (15381) . (K2T\\ and (K2^ . we obtain (KMf . 

Let us finally prove (13.251) . Let < r < f and let k = From (13.241) it 

follows that 

/ / / {k-k){x'Q + sv'o,VQ,v) dydsdCix'^^v'^) < D^A-AW^j^ll + WA- 

(5.39) 

where D4 = C2 /p dC,{xQ,VQ) and C2 is the constant that appears on the right-hand 
side of (13:241) . From (l539l) and Lemma 4.1, we obtain (K25h . □ 
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